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ABSTRACT 

ln Ecuador there is a representative number of three-story reinforced 
concrete structures with a regular configuration. It is important to indicate that many 
of these buildings have been built without adequate technical criteria, they are 
informal. Many existing buildings present high seismic vulnerability. This implies 
that these structures must be analyzed to know their behavior under the design 
earthquake established by the Ecuadorian Construction Standard (NEC). 

This article presents a methodology for the automatic generation of 
structures, based on random parameters and real information from structural 
surveys that reflect the reality of construction in our environment, in order to carry 
out a simplified seismic analysis that shows the behavior of the structures and 
based on these results take reinforcing or other actions. CEINCI-LAB, facilitates 
the development of this methodology since it allows the user to generate 
algorithms for structural analysis with different study parameters. 

From the simulation of 200 three-story structures created with random data 
and 100 real cases in Quito, with soil type D parameters, a correlation is observed 
between the drifts and the parameters of the amount of columns (area of columns 
in the first floor / area total construction), H / T (total height of the building (m) / 
Fundamental period of vibration (s)). These parameters could be used, in a 
referential way, to get an idea of the maximum drifts that can occur in a structure in 
the face of a design earthquake. 

keywords: CEINCI-LAB computer system. Reinforced Concrete frame structures. 
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GENERACION AUTOMATICA Y ANALISIS DE 
ESTRUCTURAS PORTICADAS DE HORMIGON 
ARMADO CON CEINCI-LAB 

RESUMEN 

En el Ecuador existe un numero representativo de estructuras de hormigon 
armado de tres pisos con una configuracion regular, es importante indicar que 
muchas de estas edificaciones han sido construidas sin criterio tecnico adecuado, 
es decir, que son informales y en muchos de los casos presentan alta 
vulnerabilidad ante amenaza sfsmica. Esto implica que estas estructuras deben 
ser analizadas para conocer su comportamiento bajo el sismo de diseno 
establecido por la Norma Ecuatoriana de Construccion (NEC). 

En este artfculo se presenta una metodologfa para la generacion 
automatica de estructuras, basada en parametros aleatorios e informacion real de 
levantamientos estructurales que reflejen la realidad de la construccion en nuestro 
medio, con el proposito de realizar un analisis sfsmico simplificado que muestre el 
comportamiento de las estructuras y en base a estos resultados tomar acciones de 
reforzamiento u otros. CEINCI-LAB, facilita el desarrollo de esta metodologfa ya 
que permite al usuario generar algoritmos para el analisis estructural con 
diferentes parametros de estudio. 

De la simulacion de 200 estructuras de tres pisos creadas con datos 
aleatorios y 100 casos reales en Quito, con parametros de suelo tipo D, se 
observa una correlacion entre las derivas y los parametros de cuantfa de columnas 
(area de columnas en planta baja/area total de construccion), H/T (altura total de la 
edificacion (m)/Periodo fundamental de vibracion(s)). Estos parametros podrfan 
llegar a emplearse, de forma referencial, para tener una idea de las maximas 
derivas que pueden presentarse en una estructura ante un sismo de disefio. 

Palabras Claves: Sistema de computacion CEINCI-LAB. Estructuras porticadas 
de Hormigon Armado. Generacion automatica de estructuras. 
Analisis Espectral. 

1. INTRODUCCION 

La generacion automatica de estructuras es un aspecto importante en el 
campo de la ingenierfa estructural, permite generar multiples modelos de forma 
sencilla y rapida, donde se pueden ajustar los parametros de la demanda (sismo), 
geometrfa de las estructuras (dimensiones y numeros de vanos, numero de pisos, 
alturas de entrepiso, area de construccion), dimensiones de los elementos 
estructurales (columnas, vigas, muros, etc.). 

Actualmente los programas comerciales tienen un potencial muy grande 
para el analisis estructural, pero presentan una limitacion en cuanto a modificar los 
parametros mencionados anteriormente de forma automatica, para generar 
multiples estructuras con evaluaciones de varios parametros. Sin embargo, 
algunos sistemas computacionales de codigo libre permiten realizar algoritmos 
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adaptables a las necesidades de los usuarios de forma mas flexible, en este caso 
para la generacion automatica de estructuras y el analisis modal espectral se 
emplea el sistema computacional CEINCI-LAB. 

CEINCI-LAB es una herramienta de complementacion academica para 
temas de Ingenierfa Estructural, como Analisis Matricial, Dinamica de Estructuras, 
entre otros. Adicionalmente se utiliza CEINCI-LAB para realizar investigaciones 
del comportamiento de estructuras de hormigon armado o de acero, en analisis 
simplificados de 2 dimensiones, estructuras espaciales, con disipadores de 
energia, reforzamiento de estructuras y con la aplicacion de sistemas de 
aislamiento sfsmico, entre otros temas relevantes de Ingenierfa Estructural. 

Para comprender de mejor manera la metodologfa de generacion de 
estructuras y analisis paso a paso del Metodo Espectral con el sistema de 
computacion CEINCI-LAB, el usuario puede tener como referencia los artfculos 
publicados por los primeros tres autores de este artfculo (Aguiar etal. 2019, 1 y 2). 

En Ecuador, al igual que en muchos pafses latinoamericanos, la 
construccion de edificaciones de mediana o baja altura es importante, pero 
muchas veces la informalidad genera estructuras con una alta vulnerabilidad ante 
sismos. Es por ello necesario generar metodologfas de analisis sfsmico de 
estructuras para que se pueda plantear reforzamientos de las edificaciones que no 
cumplen con los parametros mfnimos de seguridad para los ocupantes. 

En este artfculo se expone una metodologfa para el analisis sfsmico de 
edificaciones porticadas de hormigon armado, nuevas funciones de CEINCI-LAB 
para automatizar este proceso y poder analizar multiples estructuras, 
adicionalmente se estudia la relacion entre las derivas con varios parametros 
estructurales como: cuantfa de columnas y muros; periodo fundamental de 
vibracion; el parametro H/T, que corresponde a un valor resultante de la division 
de la altura total de la edificacion para el periodo fundamental de vibracion. 

2. GENERACION AUTOMATICA DE ESTRUCTURAS Y 
ANALISIS SISMICO 

Inicialmente se plantea la generacion automatica de estructuras, es decir 
que en base a ciertos parametros predefinidos se puede generar multiples 
estructuras, donde se mantiene la geometrfa de los porticos y se cambia las 
dimensiones de los elementos estructurales (vigas y columnas). Para los porticos 
en sentido longitudinal se puede calcular las matrices de rigidez y masa de la 
estructura, luego con ello se procede al analisis sfsmico mediante el Metodo 
Espectral. 

2.1. Tipologi'as de estructuras en analisis 

En algunos trabajos de investigacion se definen las tipologfas tfpicas de 
construcciones de mediana y baja altura en Quito; por ejemplo, segun la tesis de 
Cupueran en 2016, de la evaluacion estructural en 538 lotes de Turubamba en el 
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Sur de Quito, la mayorfa de las edificaciones son de hormigon armado y una 
clasificacion por numero de pisos refleja que 13% son de un piso, 41% de dos 
pisos, 34% de tres pisos y el 12% son de 4 o mas pisos. En el trabajo de Celi et al. 
en 2016, se realiza un estudio de curvas de capacidad de tipologfas estructurales; 
en este caso la informacion recopilada corresponde a 11556 edificaciones de la 
ciudad de Quito. De este trabajo se tiene clasificacion de las edificaciones en 4 
tipologfas, 16% son estructuras de un piso, 62% son parte de las tipologfas 2 y 3 
que son estructuras de 2 a 3 pisos irregulares en “L” y regulares con forma 
rectangular, 10% son estructuras de 4 a 6 pisos que tiene forma rectangular, estas 
estructuras son porticadas con vigas, columnas y losas de entrepiso de hormigon 
armado, donde se tiene mampostena de bloque. En consecuencia, se puede 
concluir de manera parcial, que, al hablar de estructuras de mediana altura en 
nuestro medio, el grupo de edificaciones de 3 pisos es representativo. En este 
artfculo se presenta inicialmente el analisis de 20 geometrfas de estructuras de 3 
pisos como se puede apreciar en la figura 1; las distancias entre ejes de columnas 
para los porticos en ambos sentidos son de 4 y 3.5 metros, lo que representa la 
tipologfa habitual de estructuras de hormigon armado en nuestro medio. 



(d) 

Figura 1 Geometrfa de las estructuras variando los vanos en sentido transversal 
con porticos tipos de sentido Longitudinal “X”, siendo estos porticos tipos 
de; a) 2 vanos; b) 3 vanos; c) 4 vanos; d) 5 vanos. 
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2.2. Modelo de analisis sismico 

De manera simplificada se realiza el analisis sfsmico de estructuras 
considerando un modelo de un grado de libertad por piso, es decir que se asume 
que el piso es rfgido y que las coordenadas principales corresponden a los grados 
de libertad horizontales, como se muestra en la figura 2. 




(b) 


Figura 2 a) Vista tridimensional de la edificacion considerando piso rfgido; b) 
Esquema simplificado de ubicacion de las masas y grado de libertad 
para el analisis sfsmico. 

2.3. Funciones de CEINCI-LAB para la generacion y analisis de 
estructuras 


Para la generacion de las estructuras se requiere el numero de vanos en 
sentido X e Y, ademas de la separacion de cada vano como ilustra la figura 3; 
para cada tipo de portico se detalla estos datos en las matrices vanosX y vanosY. 
Por ejemplo, en la matriz vanosX cada una de sus columnas corresponde a un 
portico en sentido longitudinal “X”, la primera celda de cada columna es el numero 
de vanos, luego seguido de este valor se tiene la separacion de las columnas del 
portico que define el espaciamiento de cada vano, esto para cada tipo de portico; 
de manera similar se hace en la matriz de vanosY, como se puede ver en las 
figuras 3a y 3b. 
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Figura 3 Datos de los vanos para cada estructura. (a) Sentido X. (b) Sentido Y. 
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Es importante indicar que se puede tener diferentes valores en las 
dimensiones de vanos, ademas es posible generar mas porticos de los que se 
muestra en este ejemplo. Adicionalmente para evitar problemas de definicion de 
matrices en Matlab se debe completar con ceros los espacios vacfos en las 
matrices vanosX y vanosY. Con estos datos se puede generar las diferentes 
geometrfas de las estructuras; para las estructuras en analisis se tiene 4 tipos de 
porticos en sentido “X” y 5 porticos en sentido transversal “Y”, entonces es posible 
generar las geometrfas de las 20 estructuras mostradas en la figura 1. 


Una vez asignados los vanos en sentido longitudinal de analisis “X” y en 
sentido transversal “Y”, es necesario indicar las alturas de entrepiso de la 
estructura y para ello se emplea el vector sp. En la primera parte de este artfculo, 
todas las estructuras en analisis son de 3 pisos, por lo tanto, se deben incluir 3 
datos en el vector sp. Tambien es importante conocer las dimensiones de las 
matrices vanosX y vanosY, esto se hace debido a que con esta informacion es 
posible conocer el numero total de geometrfas en analisis (en este caso 20 tipos 
de estructuras con geometrfas diversas como se muestra en la figura 1). 

Este articulo, por ser pedagogico y con el principal objetivo de ilustrar de 
manera detallada el proceso de analisis de estructuras, se presenta el codigo del 
programa en figuras; por ejemplo, la parte 1 del codigo del programa principal es 
presentado en la figura 4, con un encabezado en la parte superior izquierda, que 
presenta la secuencia del codigo. 


Codigo del programa principal. Parte 1: 


%Geometrfa del portico 
vanosX =[2 3 4 5 
4 4 4 4 

4 4 4 4 

0 4 4 4 

0 0 4 4 

0 0 0 4]; % En la primera Ifnea es #vanos y luego la separacion 

vanosY=[2 3 4 5 6 

4 4 4 4 4 

4 4 4 4 4 

0 3.5 3.5 3.5 3.5 
0 0 4 4 4 

0 0 0 4 4 

0 0 0 0 4]; % En la primera Ifnea es #vanos y luego la separacion 

sp =[2.7; 2.70; 2.70]; %lngresar un vector con la altura de pisos en (metros) 
[fil.vanosXX] = size(vanosX); 

[filY,vanosYY] = size(vanosY); 


Figura 4 Desarrollo del programa en CEINCI-LAB. 
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Las estructuras en analisis son porticadas conformadas por columnas y 
vigas de hormigon armado. La seccion de estos elementos para cada estructura 
en analisis es constante y unica, es decir una misma seccion de columna para 
todos los elementos verticales y una misma seccion de viga para elementos 
horizontales. 

Las 20 estructuras mencionadas anteriormente (figura 1), se analizan 
inicialmente con columnas de 0.1 Om x 0.1 Om y vigas bandas de 0.35m x 0.25m, 
estas columnas con una seccion tan pequena no son reales, sin embargo, esto 
representa una forma facil para generar diferentes secciones variando el lado de 
estos elementos de 5cm en 5cm, hasta llegar a una seccion de 0.55m x 0.55m (10 
tipos de columnas como se muestra en la figura 5). En consecuencia, el programa 
esta realizando el analisis de 200 estructuras. 
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Figura 5 Variacion de secciones transversales de columnas. 

Posteriormente, se debe asignar los datos del material, en este caso se 
asumen estructuras de hormigon armado con resistencia a la compresion de 210 
kg/cm 2 , para con ello determinar el modulo de elasticidad E. De manera 
conservadora, se puede asumir el valor de E = 13500* Jfc, sin embargo, es 
posible realizar una modificacion al programa para considerar un modulo de 
elasticidad diferente para cada estructura. 


Codigo del programa principal. Parte 2: 


%% Secciones % Col(varia); Vigas(35x25) (cm) 
num_tipo_col=10; 

ELEM=[ 0.1 0.1; 

0.35 0.25]; % Matriz de 2x2, En la fila 1 va base altura de columnas 
% en la fila 2 se escribe base y altura de vigas 

%% Propiedades del material 

fc =210; % La resistencia a la compresion del hormigon f'c (kg/cm2) 
E=135000*sqrt(fc); % Calcula el modulo de elasticidad del hormigon (T/m2) 


Figura 6 Datos de las secciones de vigas, columnas y materiales. 
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2.4. Datos del espectro 

Se considera el analisis en la ciudad de Quito, en base al espectro de 
aceleraciones de la NEC-15; los datos necesarios son: el tipo de suelo, la zona 
sismica, la region, el factor de amortiguamiento se asume como 5% y para el 
factor de reduccion de fuerzas sfsmicas se toma un valor de R = 5. 

De estudios de microzonificacion sfsmica (Aguiar, 2015), se menciona que 
en Quito gran parte del suelo se puede clasificar como suelos tipo D y un 
porcentaje representativo como C, siendo en menor porcentaje suelos tipo E, 
como se presenta en la figura 7. Mas adelante se requiere esta informacion para el 
analisis modal espectral. 



Figura 7 Clasificacion de suelos para Quito. (Aguiar, 2015). 

De forma didactica se asumen que las estructuras en analisis estan en 
suelos tipo D, de la ciudad de Quito, con los parametros del espectro como se 
puede apreciar en la figura 8. 


Codigo del programa principal. Parte 3: 


%% Datos para el espectro 
b=[1;1;1]; 

R=5; 

zeda=0.05; 

fip=1;fie=1; 

is=4; % Suelo D 

iz=5; % Zona sfsmica 5 = 0.4g 

ip=2; % Region Sierra 


Figura 8 Datos del espectro de la NEC-15. 

2.5. Desarrollo de analisis sismico 

El procedimiento de calculo para que sea automatico debe ser 
estandarizado, este proceso consiste en determinar la rigidez lateral de los 
porticos en sentido longitudinal “X” con cada tipo de seccion de columna, es decir 
con una determinada dimension de ese elemento como se muestra en la figura 9. 
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De forma automatica se determina el area de columnas en planta baja del portico 
en analisis y se almacena en la “matriz de celdas” denominada Acol_PB_XX. 



Figura 9 Procedimiento de calculo de los porticos tipo con diferentes secciones de 
columnas para el sentido X. 


Hasta este punto, los porticos totales calculados representan 40 matrices 
de rigidez lateral y se almacena estos calculos en la “matriz de celdas” 
denominada KL_EST, como se muestra en la figura 10 (4porticos x 10 
tipos_columnas = 40 matrices de rigidez lateral). 


KlJST _ 

;»Ui 10x4 mi _ j 

i | r .KW, | ?ft T9Ti 

2 1290*06- r»9.142- I46724U [58ST91- 

3 1878.554-. [1.1789?- 114785?- 11 7783?- 

4 [2.0330?. [2.7314?- [34289? . [4 1265?. 

5 [4*056e- (53797?- J67531?- [81266?. 

6 [7 1060?. (95322?- [1 1958?- [14384?_ 

7 [1 1702e. [15676?- [19650?- [23623e. 

8 (18205? (24358?. (30510? [36663?. 

9 [2.7054?. (3.6164«_ (45273?- [5 4383?- 
10(38705?,. [51702?- (64700?- [7 7697?- 


ICLESTO, II 
3x3 doubl? 

1 2 3 

l l 58 994 j 297018 0 4042 

2 29.7018 58.6003 -29.2821 

3 04042 -298821 288700 


Cada celda de la matriz de celdas 
KL_EST es una matriz, en este caso 
de 3x3 debido a que la estructura 
es de 3 pisos y el modelo considera 
lgdl por piso. 


Figura 10 Matrices de rigidez de los porticos tipos que componen a la estructura 
en el sentido de analisis. 
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En la figura 11 se muestra la geometrfa de los porticos en sentido de 
analisis con la numeracion de grados de libertad (para los datos de separacion de 
vanos se hace referencia la matriz vanosX). 





Figura 11 Grados de libertad en porticos de analisis. (a) Portico 1 de 2 vanos. (b) 
Portico 2 de 3 vanos. (c) Portico 3 de 4 vanos. (d) Portico 4 de 5 vanos. 


Para determinar la matriz de masas se puede asumir que la carga reactiva 
por el sismo (W) la compone la carga permanente mas la carga muerta 
correspondiente al peso de las estructuras. De forma simplificada se puede 
considerar que la dimension de la columna no sera muy relevante para el valor 
total y por ello la carga reactiva se asume como 0.8 T/m 2 para todas las 
estructuras (en futuro es posible generar variaciones de esta carga). En 
consecuencia, se requiere el area de construccion de cada piso que dependera 
unicamente de la geometrfa de las estructuras, en este caso al ser losas 
rectangulares, el area de construccion de cada piso es la multiplicacion de la 
longitud de la losa en sentido “X” por la longitud de la losa en sentido transversal, 
conforme a las geometrfas mostradas en la figura 1. Las areas de cada piso se 
almacenan en la matriz A_Const_Piso y las matrices de masa en la matriz de 
celdas denominada M_EST similar a lo mostrado en la figura 10. 

Para determinar la matriz de rigidez total de cada estructura, se multiplica 
el numero de porticos en sentido transversal por la matriz de rigidez del portico tipo 
en X, ya que los elementos y la geometrfa de estos porticos son iguales. De igual 
manera se hace para el area de columnas en planta baja. 
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Codigo del programa principal. Parte 5: 



(or I •tivaoosXX 
L4 «1 J. 

H K57( 1,111 -aJiaq (ffjau U r U |, 11,11)) ji 


Motni dr mnut para )0 c«trunura« 


Figura 12 Matrices de masas y area de cada piso de cada estructura. 

Con las matrices de rigidez y masa queda perfectamente definida cada 
estructura, con ello es posible generar el analisis modal espectral en base a la 
NEC 15 con el objetivo de calcular derivas. Es posible obtener la maxima deriva 
presente en cada estructura, indistintamente del piso en el que se genere. El 
periodo fundamental, el area de columnas en planta baja y la maxima deriva son 
almacenadas en celdas de matrices denominadas Periodos, Acol_PB_TOTALi y 
Deriva_MAX respectivamente. 


Codigo del programa principal. Parte 6: 
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Figura 13 Procedimiento del analisis dinamico. 

Para organizar esta informacion se crean las matrices Table y 
TablePeriodos que contienen informacion relevante de la estructura y de los 
resultados del analisis; por ejemplo, Table es una matriz de 200x4 debido a que 
las filas representan las 200 estructuras y las columnas de esta matriz son 
[AcolPB, AconsT, Deriva, AcolPB/AconsT] donde cada una de estas variables 
se relacionan como: 
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AcolPB, se refiere al area total de columnas en planta baja (m 2 ), siendo 
esta la suma del area de cada columna en planta baja como se muestra en la 
figura 14c; AconsT, es el area total de construccion que aporta a la carga reactiva 
del sismo (m 2 ), es decir el area sobre el nivel de planta baja como se aprecia en la 
figura 14b; Deriva, es la maxima deriva que se calcula en la estructura (%); 
Cuantia, AcolPB/AconsT, la division entre el area de columnas en PB para el area 
total de construccion (ecuacion 1). 

Cuantia de 2 A c „ ti Area de elementos resistentes a cargas laterales en PB (m 2 ) 

columnas Y,A consti Area de construccionque aporta a la cargareactiva (m 2 ) ' ' 



(a) (b) (c) 

Figura 14 (a) Vista en 3D de una estructura tipo. (b) Losas que aportan al area de 
construccion. (c) Columnas de Planta Baja que aportan al area de 
elementos estructurales. 

La matriz TablePeriodos, es una matriz de 200x2 donde las filas 
representan cada estructura, la primera columna es el periodo fundamental (s) y la 
segunda es la relacion H/T, que representa la altura total de la edificacion (m) 
divida para el periodo fundamental (s). 

El parametro H/T, es un indice de rigidez que tiene unidades de velocidad 
(m/s), esta relacion se emplea para estimar la rigidez traslacional de un edificio de 
hormigon armado, de forma general, se puede clasificar estas estructuras de 
acuerdo a la tabla 1 (Guendelman etal., 1997). 

Tabla 1 Clasificacion de estructuras en funcion del parametro H/T (Chile). 


Parametro H/T (m/s) 

Estructuras 

<20 

Muy Flexibles 

20 a 40 

Flexibles 

40 a 70 

Rigidez normal 

>70 

Rigidos 


Acorde al estudio de Moroni y Astroza (2002), en el analisis de estructuras 
despues del sismo del 3 de marzo de 1985 en Vina del Mar, se puede estimar el 
dano en funcion del parametro H/T como se presenta en la tabla 2. 
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Tabla 2 Nivel de dano esperado en funcion de H/T (Moroni y Astroza, 2002). 


Parametro H/T (m/s) 

Nivel de dano 

30 a 40 

Danos estructurales moderado 

40 a 50 

Dafio estructural leve 

50 a 70 

Dafio no estructural 

>70 

Despreciable 


Es posible observar que en la tabla 2, no se presenta estimaciones de 
dafios en estructuras de valores H/T menores a 30; no obstante, es logico suponer 
que los dafios seran mayores en los elementos estructurales y no estructurales, 
por ejemplo, Guendelman et al. (1997) indica que valores inferiores a 20 m/s 
suponen una estructura muy flexible y probablemente presentara problemas de 
cumplimiento de disposiciones de desplazamientos como se fijan en las normas de 
diseno sfsmico. Segun Cevallos et al. (2018), en estructuras con valores de H/T 
de 10 o menos se puede producir dano excesivo y hasta probablemente colapso. 

2.6. Presentacion de resultados en graficas 

En base a la informacion de resultados se puede crear graficas que 
representan el comportamiento de este grupo de estructuras. En la figura 15a es 
posible observar que la cuantfa de columnas, es decir, la relacion entre el area de 
columnas en Planta Baja respecto al area total de construccion, para estructuras 
de 3 pisos con vigas bandas, en suelos tipo D, de la ciudad de Quito, es de 0.005 
o superior para poder controlar las derivas y tener valores inferiores a 2% como lo 
establece la normativa nacional. En la figura 15b se presenta las derivas en 
funcion del parametro de altura total de la edificacion respecto al periodo 
fundamental de vibracion, en este caso valores de H/T menores a 15 generan 
derivas superiores a 2% y esto en consecuencia implicaria mucho dano ante un 
evento sfsmico. Es importante notar, que el periodo de vibracion es calculado en 
modelos numericos simplificados, sin considerar la mamposterfa; en modelos con 
mampostena el periodo disminuirfa y consecuentemente el parametro H/T podrfa 
aumentar. En el articulo de Cevallos et al. (2018), este parametro aumenta 
aproximadamente 50% al considerar la mampostena en los modelos. 

RHacion Detivas-H/T - 3 Pisos H.A. VG 0.35x0.25 



Figura 15 Graficas de resultados del analisis sfsmico (a) Derivas vs Cuantia de 
columnas. (b) Derivas vs H/T. 
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3. GENERACION DE ESTRUCTURAS A PARTIR DE BASE DE 
DATOS DE EDIFICACIONES REALES Y ANALISIS SISMICO 

En esta seccion del articulo se presenta una metodologia para generar la 
geometria de las estructuras para analisis, a diferencia del caso anterior esto se lo 
hace en base a datos reales de edificaciones; es decir, que no dependera de 
porticos tipo o de asumir dimensiones unicas de elementos estructurales. 

3.1 Configurar la base de datos de estructuras 

Para poder automatizar este proceso de generacion de estructuras, es 
necesario tener un orden especffico en la matriz de levantamiento de informacion 
de cada estructura. Por facilidad, la plantilla de toma de datos se presenta en 
Excel y esta disponible en researchgate.net/profile/Julia_Pilatasig/research junto a 
los programas necesarios para el analisis sfsmico. 

En la figura 16a se presenta el formato base del levantamiento de 
informacion de una estructura; para ejemplificar el uso de esta hoja de datos se 
presentara una estructura de 1 piso, con 3 vanos en sentido “X” y 4 vanos en 
sentido transversal “Y”. La geometrfa de la estructura se presenta en la figura 16b. 


CENTRO DE POSGRADOS ESPE 
MAESTRIA DE INVESTIGACION EN INGENIERIA CIVIL 
PROYECTO DE TITULACION 


botCni 
N° PISOS 


BOT6N2 
N° VANOS 
SENTIDO Y 



NNnutrifiN PRRA LA IXCEIEHCIA 


Datos Iniciales para Material, Cargas y Geometria. 


Porticos er 

l Sentido X 


Vigas Sentido Y 

P6rtico 1 

Pisol 

E (T/m 2 ) 

D (T/m 2 ) 

L (T/m 2 ) 

nv_x 

sv_x (m) 

np 

sp (m) 

nv_y 

sv_y (m) 

bcol (m) | 

hcol (m) 

bvg (m) 

hvg (m) 

bvg (m) 

hvg (m) 

















(a) 




(c) 

Figura 16 (a) Formato de levantamiento de informacion (HOJA_BASE). (b) Vistas 
en planta y 3D de una estructura en analisis (ejemplo). (c) Porticos “X”. 
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Para empezar, se debe hacer una copia de la HOJA BASE y colocar 
como nombre EST#, donde # corresponde al identificador de la estructura. Antes 
de llenar los datos se recomienda configurar las columnas de la hoja de datos 
siguiendo los pasos en la figura 17, en este caso dar clic en BOTON1 y colocar el 
numero de pisos (en este caso 1), siguiente dar clic en BOTON2 y escribir el 
numero de vanos en sentido transversal (sentido Y), para este caso corresponde a 
4 vanos; con esto se define el numero de columnas de la hoja en Excel para la 
geometria de la estructura. 



’i 

1 1 

T X 

UOJONl 

N'PtSOS 

hinra rf» Mai, 

•I 


(a) 




1 irtrod** 


t X 

BOT0H2 

c 4» >H(i #* i 



N* VANOS 

1 * 



STNT1DO V 






L_^_ 1 

t'iinHr 


(b) 

Figura 17 (a) BOTON1 N° PISOS. (b) BOTON2 N° VANOS SENTIDO Y. 

Los datos iniciales de material, cargas y geometria de porticos son 
registrados en la plantilla de la figura 16a; las columnas de datos son: modulo de 
elasticidad (T/m 2 ); carga muerta distribuida en la losa de cada piso (T/m 2 ); carga 
viva distribuida en la losa de cada piso (T/m 2 ); en caso de ser varios pisos las 
cargas deben ingresarse de cada piso en vertical, es decir que existira “n” filas de 
cargas en funcion de los pisos; numero de vanos en sentido “X”; longitud de cada 
vano en forma descendente (m); numero de pisos; altura de entrepiso (m); el 
numero de vanos en sentido “Y” y finalmente la longitud de cada vano en sentido 
“Y”. Estas 9 columnas son fijas e independiente del numero de pisos y vanos en 
sentido “X” o “Y”, es decir que para toda estructura se tendrfa estos datos (figura 
18a), luego es necesario las dimensiones de las columnas y vigas para poder 
definir las estructuras. En este caso, las siguientes columnas de la matriz de la 
plantilla de levantamiento de informacion tiene que ver con base, altura de 
columnas y vigas, es decir 4 columnas adicionales para cada portico en sentido 
“X”, para la estructura de la figura 16b se tiene 5 porticos en “X”, ver figura 18b. En 
este punto faltaria las dimensiones de las vigas transversales; estas dimensiones 
se colocan en las siguientes columnas, donde las dimensiones de base y altura de 
las vigas se registran para cada piso, ver los datos de la figura 18c. 


Datos Iniciales para Material, Cargas y Geometna. 

E (T/m 2 ) 

D (T/m 2 ) 

L (T/m 2 ) 

nv x 

sv x (m) 

np 

sp (m) 

nv y 

sv y (m) 

1.96E+06 

0.5 

0.1 

3 

3 

1 

2.7 

4 

4.2 

0 

0 

0 

0 

3.85 

0 

0 

0 

3.7 

0 

0 

0 

0 

2.9 

0 

0 

0 

4.35 

0 

0 

0 

0 

0 

0 

0 

0 

2.55 


(a) 
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Porticos en Sentido X 

Porticos en Sentido X 

Portico 1 

Portico 2 

bcol (m) 

hcol (m) 

bvg (m) 

hvg (m) 

bcol (m) 

hcol (m) 

bvg (m) 

hvg (m) 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.35 

0.3 

0.25 

0.35 

0.35 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0 

0 

0.35 

0.3 

0 

0 


Porticos en Sentido X 

Porticos en Sentido X 

Porticos en Sentido X 

Portico 3 

Portico 4 

Portico 5 

bcol (m) 

hcol (m) 

bvg (m) 

hvg (m) 

bcol (m) 

hcol (m) 

bvg (m) 

hvg (m) 

bcol (m) 

hcol (m) 

bvg (m) 

hvg (m) 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.35 

0.3 

0.3 

0.25 

0.3 

0.3 

0.3 

0.25 

0.35 

0.3 

0 

0 

0.35 

0.3 

0 

0 

0.3 

0.3 

0 

0 


(b) 


Vigas Sentido Y 

Pisol 

bvg (m) 

hvg (m) 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 

0.3 

0.25 


(c) 


Figura 18 (a) Datos iniciales (Primeras 9 columnas de Excel). (b) Datos de base y 
altura de columnas y vigas (Para cada portico longitudinal). (c) Base y 
altura de vigas en sentido transversal del piso 1. 

De forma general, la numeracion de columnas y vigas de un portico tipo en 
sentido “X”, se muestra en la figura 19a y para las vigas en sentido transversal, la 
numeracion de vigas se realiza similar a lo explicado en el artfculo de estructuras 
espaciales de Aguiar et al. 2019, el orden sigue como se muestra en figura 19b. 
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Figura 19 Numeracion de elementos columnas y vigas (a) Portico longitudinal, en 
sentido “X”. (b) Piso tipo de una estructura, vigas transversales. 
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Los datos en Excel, de las estructuras que se requieren analizar deben ser 
almacenados, de forma individual en archivos de texto, dentro de una carpeta 
denominada DATOS_ESTRUCTURAS. La informacion de cada estructura, 
unicamente los datos numericos (completando con ceros los espacios vacfos), 
deben ser guardadas en un archivo de bloc de notas (extension .txt) con el nombre 
de EST#.txt, donde # tiene que ver con una identificacion del numero de estructura 
(verfigura 20). 


oaioc nnurtuiAi 


Cont*e<ie 


fcST* 

l«4 


BT9 

kitto 


qp ESPE 



Figura 20 Datos de Estructuras para el ejemplo propuesto. 

3.2 Analisis sismico considerando 1 gdl por piso 


El proceso de analisis se detalla en la figura 21, los datos (EST#.txt) deben 
ser almacenados en una carpeta denominada DATOS_ESTRUCTURAS. 
Posteriormente se debe ejecutar el programa de CEINCI-LAB 
Modelos_JEP_Automatizacion_1GDL.m que genera la estructura y realiza el 
analisis de cada modelo, de manera similar a lo que se indico en la seccion 2 del 
presente artfculo; en este caso la diferencia es que la geometrfa, cargas, 
materiales, dimensiones de columnas y vigas se encuentran previamente definidas 
en los archivos EST#.txt. Al ejecutar este modelo se genera automaticamente una 
carpeta RESULTADOS_ESTRUCTURAS_1GDL con los resultados de los analisis 
sfsmicos realizados, en esta carpeta se generan archivos de textos que 
almacenan los resultados de cada estructura, que son denominados 
RES_1GDL_EST_#.txt. Finalmente, para poder trabajar con los resultados se 
debe ejecutar el programa Resultados_JEP_Automatizacion_1GDL.m que 
genera graficas como se aprecia en la figura 23, que son similares a las 
presentadas en la figura 15, pero con los datos de cada estructura. 



G«om«tna, 

Ntatcriaics. 

Srcoonrv. 

Car*»s, 

Etc. 



# ModekJSjtf _Autcm«nrioor JGLH 

for I = 1 ; nom_£itructur»s • 

l£dl por 
pko 


Or?nv?»v, Ppriodo fundam«nt»i 
(T), Aitura ndifmo |H), 
pararrvrtro HA, *tc. 

* end 


4 itMUtadMjEF.Auiomataicicn.lGOl 



4 

W.. 


Retuiudos 
de cada 
estrudura 


Figura 21 Proceso de analisis sfsmico de las estructuras. 
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Para el calculo de la matriz de rigidez, se considera un modelo similar al 
presentando en la figura 2, es decir, un modelo con 1 grado de libertad por piso; 
se realiza un analisis en sentido X y posteriormente se cambian los grados de 
libertad para hacer un analisis en sentido transversal Y (modelos de la figura 22). 




Figura 22 Modelos de 1 gdl para analisis sfsmico (a) Analisis en sentido 
longitudinal “X”. (b) Analisis en sentido transversal “Y”. 

La codificacion del programa empleado en esta seccion es similar al 
presentado anteriormente; en la figura 23 se puede apreciar parte del codigo del 
programa Modelos_JEP_Automatizacion_1GDL.m que difiere de lo presentado 
en otros artfculos sobre el analisis sismico con CEINCI-LAB. En este caso se 
puede notar que al inicio se debe dar el nombre de la carpeta 
RESULTADOS_ESTRUCTURAS_1GDL que almacenara los resultados del 
analisis de cada estructura, seguido del numero de estructuras que estan en la 
carpeta de datos, para este caso 50 estructuras con datos reales correspondientes 
a viviendas de Hormigon Armado de 1 a 3 pisos y una estructura de 4 pisos en la 
ciudad de Quito. 

El analisis sfsmico de este grupo de estructuras es secuencial, es decir 
desde la estructura 1 a la nEst, es por ello que se presenta un bucle for ... end en 
el cual se alberga el codigo que realiza los calculos del analisis. En la figura 23 es 
posible notar que la informacion de cada estructura se almacena en una matriz 
denominada EST, la dimension de esta matriz varia en funcion de la geometria de 
cada estructura (numero de vanos y numero de pisos), tambien varfan las 
dimensiones de las matrices que se involucran en el analisis, por ejemplo las 
matrices de rigidez y de masa depende del numero de pisos de cada estructura, 
por ello antes de finalizar el bucle se coloca un comando de Matlab “clearvars 
except ij” que elimina todas las variables del espacio de trabajo (rigidez, masa, 
np, sv, etc.) con excepcion del identificador de la estructura en analisis 
denominado ij, esto se lo realiza para evitar problemas de dimensiones de 
matrices en los calculos y tener la seguridad de que cada analisis es realizado con 
los datos que corresponden a cada estructura. 

Las columnas de la matriz de datos estan ordenadas de forma que se 
pueda aplicar un programa para automatizacion del analisis, una nueva funcion de 
CEINCI-LAB es Secciones_X_Y_EST con la cual se obtienen las dimensiones de 
las columnas y vigas de cada portico (base y altura de los elementos, asumiendo 
elementos rectangulares) para el analisis en sentido X y Y. 
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Modelos JEP Automatizacion 1 GDL.m: 

%% Carpeta en la cual se van a grabar los archivos de resultados 
mkdir RESULTADOS_ESTRUCTURAS_1 GDL 

nEst = 50; % Numero de estructuras en analisis (Mismo # de archivos txt) 

for ij=1 :nEst 

%% Leer cada archivo con los datos de cada estructura 
EST = load(['DATOS_ESTRUCTURAS\EST' num2str(ij) ’.txt']); 

E=EST(1,1); % Modulo de elasticidad, T/m2 

Cmuerta=nonzeros(EST(:,2)) ; % Carga muerta en cada piso, T/m2 

Cviva=nonzeros(EST(:,3)) ; % Carga viva en cada piso, T/m2 

nvx=EST(1,4); % Numero de vanos en X 

svx=EST(1 :nvx,5); % Longitud de cada vano en X, m 

np=EST(1,6); % Numero de pisos 

sp=EST(1 :np,7); % Altura de entrepiso, m 

nvy=EST(1,8); % Numero de vanos en Y 

svy=EST(1 :nvy,9); % Longitud de cada vano en Y, m 

% Numero de elementos de un portico en sentido X 

numColX=np*(nvx+1); 

numVigasX=np*(nvx)*2; 

% Numero de elementos de un portico en sentido Y 

numColY=np*(nvy+1); 

numVigasY=np*(nvy)*2; 

% Secciones de cada elemento de cada portico 

[SeccionPorticosX,SeccionPorticosY]=Secciones_X_Y_EST(np,nvx,nvy,EST); 
Analisis sismico en sentido X considerando 1GDL por piso 

Metodologia similar a la presentada en la seccion anterior 
Analisis sismico en sentido Y considerando 1GDL por piso 

Metodologia similar a la presentada en la seccion anterior 
ResultadosGeneral=[np, Acol_PB_TOTAL, A_Const_TOTAL, HTedificio, 
maxDervX, maxDervY]; 

[f 1 ,c1 ]=size(ResultadosGeneral); 

ResultadosX=[Masa_Traslacional, KL_T_X, TX, qteX]; 
[f2,c2]=size(ResultadosX); 

ResultadosY=[KL_T_Y, TY, qteY]; 

[f3,c3]=size(ResultadosY); 

ResultadosDESP=[qtiX', DervX', qtiY', DervY']; 

[f4,c4]=size(ResultadosDESP); 
filas=max([f1 ,f2,f3,f4]); 
columnas=c1 +c2+c3+c4; 

Resultados=zeros(filas,columnas); 

Resultados(1 :f1,1 :c1)=[ResultadosGeneral]; 

Resultados(1 :f2,c1 +1 :c1 +c2)=[ResultadosX]; 

Resultados(1 :f3,c1 +c2+1 :c1 +c2+c3)=[ResultadosY]; 

Resultados(1 :f4,c1 +c2+c3+1 :c1 +c2+c3+c4)=[ResultadosDESP]; 

%% Grabar los archivos 

save(['RESULTADOS_ESTRUCTURAS_1 GDL\RES_1 GDL_EST_' num2str(ij) 
'.txt'],['Resultados'], -ascii'); %.DATO: Direccion carpeta 
clearvars -except ij 

end 


Figura 23 Fragmento del programa principal para el analisis sismico. 
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3.3 Presentacion de resultados y generacion de graficos 

En los archivos de resultados de cada estructura RES_1GDL_EST_#.txt, 
que se almacenan en la carpeta RESULTADOS_ESTRUCTURAS_1GDL, se 
presenta la informacion para cada estructura, ordenada de la siguiente manera: 
primera columna es el numero de pisos de la estructura; el area de la seccion 
transversal de todas las columnas en planta baja (m 2 ); area total de construccion 
(m 2 ); altura total de la edificacion (m); maxima deriva presentada en el analisis en 
“X” (%); maxima deriva presentada en el analisis en “Y” (%); matriz de masa 
(T*seg 2 /m); matriz de rigidez total en sentido “X” (Tm), periodos de vibracion del 
analisis en sentido “X” (s); desplazamientos elasticos en cada piso del analisis en 
“X” (m); de manera similar en sentido “Y” se presenta matriz de rigidez, periodos 
de vibracion y desplazamientos elasticos del analisis; continuan los resultados con 
desplazamientos inelasticos (m) y derivas (%) en cada piso para el analisis en 
sentido “X”, finalmente estos mismos parametros para el analisis en sentido “Y”. 

El programa Resultados_JEP_Automatizacion_1GDL.m genera las 
graficas con los archivos de resultados, su metodologia consiste en extraer los 
datos resultantes del analisis sismico e identificar los parametros del archivo 
(matrices y vectores), en funcion del numero de pisos de cada estructura y 
posteriormente realizar las graficas que el usuario quiera para hacer una discusion 
de analisis, la parte principal del codigo del programa se muestra en la figura 24. 


Resultados JEP Automatizacion IGDL.m: 


nEst = 50; % Numero de estructuras en analisis (Mismo # de archivos txt) 

for ij=1 :nEst 

%% Leer cada archivo con los datos de cada estructura 
EST{ij,1} = load(['RESULTADOS_ESTRUCTURAS_1 GDL\RES_1 GDL_EST 
num2str(ij) '.txt']); %.DATO:Direccion carpeta 
np(ij,1)=EST{ij,1}(1,1); 

Acol_PB_TOTAL(ij,1 )=EST{ij,1 }(1,2); 

A_Const_TOTAL(ij,1 )=EST{ij,1 }(1,3); 

HTedificio(ij,1)=EST{ij,1}(1,4); 
maxDervX(ij,1 )=EST{ij,1 }(1,5); 
maxDervY(ij,1)=EST{ij,1}(1,6); 

TX(ij,1 )=EST{ij,1 }(1,6+2*np(ij,1 )+1); 

TY (ij, 1 )=EST{ij,1 }(1,6+1 +3*np(ij,1 )+1+1); 

% Calculos 

Cuantia_Col(ij,1)=Acol_PB_TOTAL(ij,1)/A_Const_TOTAL(ij,1); 

HTX(ij, 1 )=HTedificio(ij,1 )/TX(ij, 1); 

HTY(ij, 1 )=HTedificio(ij, 1 )/TY(ij, 1); 

end 


Figura 24 Fragmento del programa principal para resultados. 

Se puede notar en la figura 25, que los resultados de analisis de casos 
reales son similares a los presentados en la figura 15, correspondientes a 
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estructuras generadas aleatoriamente. No obstante, la cuantia rmnima de 
columnas en este caso deberia ser 0.004 para cumplir con las derivas maximas 
permitidas de la NEC-15, esto se debe a que las estructuras reales representan 
modelos de estructuras con vigas bandas y vigas peraltadas. 


Cuantia de columnas para suelo 0 en Quito 
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Figura 25 Graficas de resultados del analisis sfsmico (a) Derivas vs Cuantfa de 
columnas. (b) Derivas vs H/T. 

3.4 Analisis complementarios sobre las cuantias de columnas 


Se realiza el analisis de estructuras de 3 pisos de Hormigon Armado 
ubicadas en suelo tipo D de la ciudad de Quito, considerando vigas bandas de 
35x25, vigas peraltadas de 30x40 y otras estructuras con vigas de 25x35, con la 
metodologfa de generacion aleatoria de estructuras en base a parametros 
predefinidos (seccion 2). En la figura 26, se presenta los resultados mas 
relevantes de estos analisis junto a los resultados de los analisis de estructuras 
reales (seccion 3), es posible observar que las estructuras aleatorias engloban a la 
tipologfa de estructuras reales y hacen posible visualizar una tendencia 
exponencial decreciente entre las variables de maximas derivas (%) y cuantfa de 
columnas, destacando como valor referencial una cuantfa mmirna de 0.004 para 
obtener derivas inferiores a 2%. 
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Figura 26 Derivas vs Cuantfa de columnas. 
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Es posible generar analisis de este parametro en funcion del tipo de suelo, 
sin embargo, al ser estructuras de 3 pisos donde la mayorfa de estas estructuras 
tienen periodos hasta 0.6 segundos aproximadamente, entonces el tipo de suelo 
no genera un cambio representativo, por ejemplo, en la figura 27a se puede 
observar una tendencia similar a la presentada en la figura 26, siendo adecuado 
una cuantfa de 0.004 para controlar las derivas. En la figura 27b se muestra como 
influye el tipo de suelo en las estructuras con geometrfa real en la ciudad de Quito 
(seccion 3 de este artfculo), se aprecia que existe poca variacion entre los 
resultados para suelo C y D, esto se debe a que para el rango de periodos de 
estas estructuras de 1 a 3 pisos los espectros de la NEC 15 para suelos C y D 
presentan una forma similar (magnitudes de aceleraciones), las estructuras en 
suelos tipo E varian debido a que las aceleraciones del espectro NEC 15 para 
suelo E es menor en este rango de periodos. 



Cijarba de cofcmnas (•) Cuarcta da cotumno* (•) 

(a) (b) 

Figura 27 Derivas vs Cuantfa de columnas. 

4. CONCLUSIONES 


En este articulo se presenta una metodologfa para generar estructuras de 
manera automatica, con la combinacion de parametros predefinidos, y tambien 
con informacion de una base de datos de estructuras pertenecientes a Quito. 
Esto permite optimizar el tiempo de calculo y generar nuevas investigaciones en 
las cuales se estudie el impacto de los parametros estructurales y de la demanda 
sfsmica, en el comportamiento de las edificaciones. Los resultados presentados 
son parciales y pertenecen a la tesis de Maestrfa de Julia Pilatasig. 

El sistema de computacion CEINCI-LAB permite realizar modificaciones a 
los codigos presentados, para adaptarse a las necesidades de cada usuario; 
debido a esta flexibilidad y a que es de libre acceso representa una ventaja para 
estudiantes, profesionales e investigadores en el area de estructuras. 

Del analisis sfsmico de 200 estructuras de hormigon armado de 3 pisos, 
con vigas bandas (35x25), 400 estructuras con vigas peraltadas (25x35 y 30x40) 
y de 100 casos adicionales de estructuras reales de 1 a 3 pisos con diferentes 
geometrias y dimensiones de elementos estructurales, simuladas sobre suelo tipo 
“D” en la ciudad de Quito, se puede inferir que existe una correlacion entre las 
derivas de piso y la cuantfa de columnas, para cumplir con derivas inferiores al 
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2%, que es el Ifmite de la NEC-15, se recomienda una cuantia de columnas 
minima, es decir, una relacion de area de la seccion transversal de todas las 
columnas en planta baja respecto al area total de construccion superior a 0.004. 
Este valor es referencial y solo es aplicable en base a las condiciones estudiadas, 
mas adelante se presentara resultados con mayores simulaciones. Respecto al 
parametro H/T, con valores menores a 15 m/s las estructuras podrfan ser 
categorizadas como muy flexibles, con derivas superiores a 2%, por lo tanto, 
podrfan presentar mucho dano ante un evento sfsmico. 

Los parametros presentados en este artfculo, para predecir el 
comportamiento sfsmico, son faciles de determinar, de manera que pueden 
representar un buen valor referencial para conocer de manera simplificada las 
derivas de las estructuras ante un sismo de diseno e incluso se podrfan emplear 
para evaluaciones rapidas de vulnerabilidad estructural, luego en base a estos 
resultados se puede considerar estudios mas detallados en las estructuras que 
considere el ingeniero Estructural. 
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